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Abstract 

Depinning of an interface from a random self-affine substrate with rough- 
ness exponent Cs is studied in systems with short-range interactions. In 2D 
transfer matrix results show that for C^s < 1/2 depinning falls in the univer- 
sality class of the flat case. When ()s exceeds the roughness (Co = 1/2) of the 
interface in the bulk, geometrical disorder becomes relevant and, moreover, 
depinning becomes discontinuous . The same unexpected scenario, and a pre- 
cise location of the associated tricritical point, are obtained for a simplified 
hierarchical model. It is inferred that, in 3D, with Co = 0, depinning turns 
first-order already for Cs > 0. Thus critical wetting may be impossible to 
observe on rough substrates. 

PACS numbers: 64.60.Ak, 68.45.Gd, 36.20.Ey, 75.60.Ch 
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Wetting and depinning phenomena occnr when the interface between two coexisting 
phases nnbinds from an attractive snbstrate [Q]. In 2D, where an Ising interface can be well 
represented by a directed path [^, critical wetting occnrs as a rule, and hrst-order wetting is 
predicted only for special setups |^. In 3D, on the other hand, numerical simulations have 
shown critical, hrst-order and tricritical wetting for the Ising model. The same phenomena 
in random media have attracted a lot of attention recently. Many works M considered the 
effects of quenched disorder due to bulk impurities on interface depinning. Other studies ad- 
dressed disorder restricted to a geometrically smooth surface (chemical surface disorder) [Q . 
Both bulk and chemical surface disorder may modify the universality class of the unbinding 
transition, but generally not its continuous, second-order character in 2D. JI| 

In the present Letter we address the role of geometrical disorder due to wall roughness in 
determining the nature of the wetting transition. Disordered geometry is most amenable to 
theoretical treatment when characterized by simple scaling laws, like in the case of self-affine 
or fractal substrates. So far, the effects of such roughness were seldom discussed, mostly 
in connection with complete wetting [^|-p. In spite of this, real substrates with self-affine 
geometry are met in many situations and have been the object of recent experiments [Q. 
For these substrates the average transverse width of the wall scales with the longitudinal 
length X as X**® (0 < ^5 < 1). 

We show that in the 2D Ising model the nature of the depinning transition changes 
drastically upon increasing the self-affine roughness exponent C,s- For (C5 ^ 1/2 depinning 
remains continuous as in the flat case, and disorder is irrelevant. For > 1/2 an unusual 
discontinuous depinning occurs. Since the intrinsic roughness of a 2D interface is (/o = 1/2, 
it is natural to expect Cs = 1/2 as precise threshold. The emerging scenario is that of geo- 
metrical disorder determining a tricritical phenomenon for depinning. We are able to locate 
the tricritical point by accurate renormalization group (RG) calculations on a hierarchical 
model, which further suggest = 1/2 as the tricritical threshold in 2D. Extrapolation of 
our results to the experimentally most relevant case of 3D, where </o = 0 in ordered bulk, 
suggests that a minute disorder in substrate geometry could suppress critical wetting in 
favor of a hrst-order transition. 




FIG.l. Example of rough substrate boundary (continuous) and interface (dotted) conhgu- 
r at ion. 

Let X and y be integer coordinates of points on a square lattice. A “waii” is given by 
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a self-avoiding path directed according to the positive x axis (Fig. 1). To each wall step 
parallel to the x axis, a height Sx, eqnal to the ordinate of its left end, x, is associated. For 
simplicity we consider only wall conhgurations obeying the restriction: Sx = Sx — Sx+i = ±1. 
In order to generate wall geometries with a preassigned ronghness exponent (s, we nsed a 
randomized version of an algorithm dne to Mandelbrot m- Given a wall prohle, the interface 
can assnme conhgurations obeying hx — hx+i = 0, ±1 and hx > Sx, if hx is the height of the 
horizontal step at x (Fig. 1). The interface is like a (partially) directed walk or polymer, 
and its Hamiltonian is 



Ff 'y ( [^ ( f T \zx Zx-\-\ T (f) 

X 

where Zx = hx — Sx, and e,u> 0. According to (1), at temperature T, fugacities a; = 
and k = are associated to each (horizontal or vertical) step of the walk, and to each 
horizontal step on the wall, respectively. With a wall prohle covering a distance X the 
partition function is 

Zx= E (2) 

all walks 

where the sum is restricted, e.g., to walks from the origin (0, 0) to any point {X, y), with y > 
Sx- L and I indicate total length and number of horizontal steps on the wall, respectively. 
Zx is a functional of the wall prohle. In order to calculate it we use transfer matrices: 

(Tsj,)m,n ^\hn,mr-Sx~h ’ (3) 



where m and n range on the allowed Zx and Zx+i, respectively. The partition function thus 
becomes 






Efll'n.) Up) 

l,p \a:=0 / l^p 



( 4 ) 



where, with the left end of the interface grafted at the origin, 0o(p) = ^p,o- A wall prohle 
corresponds to a sequence of factors Ti, T_i in the product of eq.(4). Zx is thus given 
asymptotically in terms of the largest Lyapunov eigenvalue H 



^max 



= lim 

X^oo 




X 



( 5 ) 



We verihed that diherent long T sequences, i.e. wall prohles, lead to the same eigen- 
value within good accuracy. Thus, the quenched free energy is liaXmax{^,k,Cs) = 
limx^oo where the bar indicates quenched averaging over wall prohles. In the 

random context depinning is most efficiently detected by studying the behavior of quantities 
which can be directly related to the components of (j)x- Examples are the average probability, 
Pq{uj, k, (s) = limx^oo y ^x=o 0i(O)/I 10*1 P; that the horizontal step lies on the wall, and the 

average distance of the interface from the substrate, (z) = lim^^oo x '^xZo Sz ^0s(^)/l 10*1 P 
m]. For our determinations we used up to 50 independent prohles with x < for which 
the components of (px were computed up to a distance from the wall Zmax — 2 ■ 10^. 
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TABLE I. Values of V’ and kc sX ui = 1/2 from fits of {z). 



Cs 


0 


kc 


0 


1 


4/3 “ 


1/3 


1.01 ± 0.01 


1.772 ± 0.001 


2/5 


1.01 ± 0.02 


1.828 ± 0.001 


1/2 


1.03 ± 0.04 


1.908 ± 0.003 



“ Exact results for a flat substrate [1]. 



Rather than considering variations of Pq or {z) along curves parametrized by tempera- 
ture, we choose to follow a; = const, lines. In the case (s ~ 1/2, e.g., upon approaching 
k = kc from above with 0 < a; < 1, (z) is well htted by (z) = A{k — kc)~^^ + B , with -0 always 
compatible with the exactly known flat wall value of 1 [^. kc of course depends on oj and (s- 
Some "0 and kc determinations are reported in Table I for u> = 1/2. For (s ~ 1/2 disorder 
in the wall geometry does not appear to lead to a new universality class for depinning. By 
treating a different model with continuum many-body techniques, Li and Kardar found 
second-order depinning with 0 = 1 for all < 1. At variance with this conclusion we 
hnd here that the situation drastically changes for (s ~ 1/2. In this range (z) has a much 
steeper, abrupt rise at kc., so that the previous ht becomes clearly inappropriate. 

Further strong evidence of a change of the nature of the transition at </s — 1/2 comes 
from the behavior of Pq (Fig. 2). While for </s ~ 1/2, Pq ~ (k — kc)^, for k — > with p — 1, 

as with flat substrate |^, for ^5 > 1/2 a discontinuity in Pq shows up (Fig. 2), becoming 
sharper and sharper as hnite-size effects are reduced. The amplitude of the discontinuous 
jump in Pq also increases with (/s- 




_ k 

FIG. 2. The probability Pq as a function of /c at a; = 1/2 and for different C,s values. 
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Thus, for ^5 > 1/2 substrate roughness is relevant and, moreover, drives the depinning 
transition hrst-order. (^5 = 1/2 is the natural candidate as border value between continuous 
and discontinuous regimes. Indeed, for > 1/2 the wall roughness exceeds the roughness 
Co = 1/2 P of the interface. 

In 2D, hrst-order depinning is quite unexpected in the context of interfacial phenomena. 
Only two special ways of obtaining it have been conceived so far, by introducing either 
an attractive defect line in the bulk []^, or longitudinally fully correlated disorder []^ 



Here hrst order is caused by sufficiently strong geometrical surface disorder, which also 
reveals opposite in its effects to its chemical counterpart. Indeed, while higher roughness 
induces hrst-order, in the defect line case chemical surface disorder drives the transition 
back continuous [Q. This latter ehect is certainly what one would expect at hrst sight 

Like in the special examples of refs. 
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on the basis of experience with phase transitions 
|T^,[T^, our hrst order depinning needs not be accompanied by oh-coexistence prewetting 
phenomena. 

Extrapolation to 3D of our hndings is natural and has remarkable and unexpected conse- 
quences. Systems with a dominance of short-range interactions are, e.g., metallic substrate- 
adsorbates or, even more, type-I superconductors m- Since typically interfaces in pure 
systems are only logarithmically rough in 3D ((/o = 0), a minimum of substrate roughness 
should be sufficient to give hrst-order wetting. This ohers a further possible explanation 
of the fact that critical wetting is so elusive from the experimental point of view. A most 
recent work predicts critical wetting for superconductor interfaces [jT^. For such interfaces 
Co is not known, unfortunately, but could be rather small, if not zero. This means that 
special care in using smooth substrates should be exerced, in order to observe the predicted 
phenomenon. 





FIG. 3. Construction rule of the DHL (level n = 0 to level n = 1) and schematic picture 
of the lattice at level n, with the four n — 1 level units. A wall conhguration (heavy) 
crossing the left units, and two polymer conhgurations (dotted) are reported. With such 
wall conhguration eq. (6) applies. 



By reinterpreting to and k as monomer fugacity and Boltzmann factor for contacts, 
respectively, our model describes polymer adsorption [j^. Apart from a change in the 
ensemble {Z = Z^x^v), the transfer matrices are the same. Criticality {Z dominated by 
inhnite length polymer conhgurations) implies Xmax = 1- We hnd that Amaa:(l/2, fc, Cs) = 
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1 for all k < kc{l/2Xs), where kc was defined above. For k > kc criticality occurs at 
oj < 1/2, indicating that the polymer is adsorbed Indeed a; = 1/2 marks criticality 

for the polymer in the bulk. The dependence of k^. on C,s shows that adsorption becomes 
more difficult with increasing roughness. One can also show that the fraction of monomers 
adsorbed on the wall should have the same singular behavior as Pq, when moving at a; = 1/2. 
Thus, like interface depinning, polymer adsorption undergoes a change from second to hrst 
order upon increasing X- 

To gain more insight into this change, we stick to polymer language and consider a 
simplihed model of adsorption dehned on diamond hierarchical lattice (DHL) (Fig. 3). Self- 
avoiding paths on DHL have often been used to mimic directed polymers in 2D |T^. A wall 
joining the two ends is obtained as follows: at level n = 0 of DHL construction the wall 
always coincides with the unique existing bond. For n = Umax > 0, the wall is determined 
by backward iteration. Starting from rimax, at each level, n, we choose whether the wall 
passes through the left or right DHL units of level n — 1, with probabilities 1 — A and A, 
respectively, and so on. If, e.g., we put A = 0, the process is deterministic and we create 
a single wall coinciding with the left border of the lattice. For A = 1/2 we generate with 
equal probability all possible walls through the lattice. Given a wall, we consider a polymer, 
with partition Zn, joining the ends of the lattice and laying, e.g., to the right of the wall. 
As in eq. (2), oj and k are step and wall contact fugacities, respectively. For A = 0 we deal 
with a polymer attracted {k > 1) by the left DHL border, which plays the role of a “flat”, 
deterministic substrate. When A rises, the “roughness” of the now random wall increases. 
Transverse hills and valleys are felt more and more by the polymer. For a given wall, one 
can compute the polymer partition function iteratively, using Zh^n+i = and 



2^n+l — Zn^\Zn^2 + 



( 6 ) 



or 

Zn+l = Zn^\Zn^2 ( 7 ) 

Eq.(6) or (7) is chosen, according to whether, at level n + 1, the n-th level diamonds, 1 
and 2, crossed by the wall, are the left or the right ones, respectively. is clearly the 
“bulk” partition function on DHL at n-th level, in absence of wall. Initial conditions are 
Zq = kui and Zi,^ = oj. Considering hrst A = 0, eq.(6) induces a two parameters RG 
mapping by putting 2^1 = uj' and Zi = uj’k’ . Clearly the transformation of Zi^^n implies 
that the bulk criticality condition is a; = 1/2 as in the Euclidean case. The value k = 1 
is marginally unstable for A; > 1 and separates adsorbed {k > 1) from desorbed {k < 1) 
regimes. Thus, for A = 0 a positive attraction is always sufficient to adsorb the polymer. 
Eq. (6) applies for all n and the problem has a nontrivial hxed point with Z^ = 1/2 and 
Z* = 1/2. Po{oj = 1/2, fc, A = 0) can be extrapolated by iteration. Due to marginality, 
Pq starts rising with zero slope, but continuously, for k = 1 (Fig. 4). For A > 0 the 
ZnS become random variables and we must iterate their probability distribution, V^. This 
cannot be done exactly. However, starting from Vq = 6{Z — ku>), we follow the distribution 
at level n by iteratively sampling it. From a large number (~ 10®) of Z values distributed 
according to Vn-i) we generate a sample distributed according to Vn by choosing many pairs 



6 



of Z values and obtaining from each pair a new value of Z according to eq.(6) or (7), with 
probabilities 1 — A and A, respectively. The resulting 2’s constitute a sample of Vn- This 
procedure could be iterated for n < 40 with extremely high accuracy. 




_ k 

FIG. 4. Pq on the DHL at cn = 1/2 for various A values. The dashed lines mark the 
discontinuities for A > 1/2. 

Pq at cn = 1/2 is plotted as a function of k in Fig. 4 for different A’s. For 0 < A < 1/2, 
Po rises continuously from zero as c{k — kc). Thus, p = 1 for 0 < A < 1/2. The slope 
c and kc are both increasing with A. Apart from c = 0 at A = 0, due to the accidental 
marginality, such behaviour reproduces what is observed in 2D when < 1/2. For A < 1/2 
the polymer is more “rough” than the wall. Only when A = 1/2 the latter has the same 
freedom to develop through the DHL as a polymer has within the “bulk”. On the other 
hand, A > 1/2 corresponds qualitatively to > 1/2, because the polymer feels more and 
more the wall limiting its options when developing through the DHL. For A = 1/2 we have 
evidence that c = oo, with a still continuous transition. This inhnite slope indeed anticipates 
a sharp discontinuity in Pq for A > 1/2. So, the hierarchical model contains ingredients 
reproducing, at least qualitatively, the scenario emerging for the Euclidean model, and gives 
a suggestive indication of the way in which continuous transitions switch to hrst order at 
the expected threshold (s = 1/2. 

The dependence of kc on A mimics that on (s in 2P, and further motivates the cor- 
respondence between A and (s in the two cases. We stress that, since all paths have the 
same length, and there is no natural recipe for dehning a transversal distance on DHL, the 
notion of roughness must always be mediated in some way: here we can link roughness to 
A. Our hierarchical model provides a remarkable example of the tricritical transition we are 
dealing with in this paper, and allows an essentially exact determination of its location and 
properties. Moreover, the threshold we hnd at A = 1/2 is strongly suggestive of the precise 
location of the tricritical point in the Euclidean case. 

Summarizing, in 2D, interface depinning or directed polymer adsorption on rough sub- 
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strate with (s < 1/2 are continuous and in the same universality as in the flat case. For 
Cs > 1/2 roughness is relevant and, moreover, the transitions acquire an unusual, discontin- 
uous nature. This result, not anticipated so far |0], warns that some continuum approaches 
may not be able to catch the correct physics of depinning from rough substrates. We expect 
similar properties, and the same threshold (s = 1/2, for directed polymer adsorption on a 
self-affine surface in 3D. Indeed, directed polymers have no upper critical dimension, and 
for them = 1/2 in all d. Although the single polymer adsorption regime is not easily ac- 
cessible experimentally, we believe that our results should be relevant for stretched polymers 



El. 

In 3D the interface of a pure system typically has Co = 0 |^. We thus conjecture that 
depinning occurs discontinuously as soon as Cs > 0, in cases when the interactions are pre- 
dominantly short-range, like in metallic systems. Our results also bear on the observability 
of the recently predicted critical wetting in the case of type-I superconductors |T^, which 
is perhaps the most strict physical example of short-range interface-substrate interactions. 

Numerical calculations have been partly supported by CNR within the CRAY project of 
Statistical Mechanics. We thank Joseph Indekeu for valuable criticism and suggestions and 
Mehran Kardar for stimulating discussions. 
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